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Abstract 

The atomic and band structures of the (6, 0) zigzag carbon nanotube at its axial elongation are calculated by semiempirical molecular 
orbital and by tight-binding methods. The ground state of the nanotube is found to have a Kekule structure with four types of bonds 
and difference between lengths of long and short bonds of about 0.005 nm. The structural phase transition is revealed at w 9% 
elongation, resulting in a quinoid structure with two types of bonds. This structural phase transition is followed by the transition 
from a narrow gap to moderate gap semiconductor. Validity of the semiempirical PM3 method is discussed. 

Keywords: Zigzag carbon nanotube, Molecular orbital calculations, Kekule structure, Phase transition, Tight-binding 
approximation, Band gap 
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1. Introduction 

The studies of structural, electronic and elastic properties of 
carbon nanotubes (CNTs) are actual in connection with per- 
spectives of their applications in nanoelectronic devices and in 
composite materials [ 1 ] . These properties are also of fundamen- 
tal interest, particularly for physics of phase transitions. For ex- 
i ample, such structural phase transitions as the commensurate- 
incommensurate phase transition in doublewalled CNTs |2j |3J 
and the spontaneous symmetry breaking with formation of cor- 
rugations along nanotube axis [4] have been considered. 

The possibility of Peierls transition in CNTs was first con- 
sidered in [5]. As a result of this transition, metallic CNTs 
become semiconducting at low temperature and Peierls distor- 
tions of the nanotube lattice lead to a Kekule structure (see 
Fig. la). The Kekule structure can arise in a (m, n) CNT in 
the case that m — n is multiples of 3 J6] . The Peierls gap 16- 
8], Peierls transition temperature 

iSiiHIl 

and bond lengths 

changes JfJ Q, due to Peierls distortions were estimated for 
(n,n) armchair CNTs. Static twist deformation of the (5,5) 
armchair nanotube lattice at the Peierls transition was also pre- 
dicted 111, 12h . The other type of Peierls distortions was found 
for the semiconductor chiral CNTs [3 ■ The Kekule struc- 
ture of the ground state was found for the finite length armchair 
CNTs using density functional theory 1 15 - 17| and semiempir- 
ical molecular orbital 1U7I1 calculations. X-ray crystallographic 
analysis of chemically synthesized short (5, 5) CNTs shows the 
Kekule bond length alternation pattern for their structure 01711 . 
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Although the Kekule structure can arise in («, 0) zigzag 
CNTs, if n is multiples of 3, such possibility has not yet been 
considered. The density functional theory calculations per- 
formed for short zigzag CNTs show that pattern of short and 
long bonds arise II 1 8L Il9ll . however this pattern cannot be as- 
signed to the Kekule structure (as it is defined in [6, 20]). Previ- 
ously 11 1 311 we have revealed the Kekule structure of the ground 
state of (5,5) armchair CNT by semiempirical molecular or- 
bital calculations and found two deformational structural phase 
transitions connected with spontaneous symmetry breaking and 
controlled by uniaxial deformation of the CNT. 

In the present Letter we consider the possibility of the Kekule 
structure of the ground state and structural phase transitions 
in (3n,0) zigzag CNTs on the example of the infinite (6, 0) 
CNT. We present the semiempirical molecular orbital PM3 02 ill 
calculations of the internal energy and atomic structure of a 
uniaxially deformed infinite (6,0) CNT up to 12% elonga- 
tion. For calculated atomic structure we use simple tight- 
binding (Hiickel) 02211 calculations to gain qualitative insight 
into changes of the electron energy band structure of a uniaxi- 
ally deformed (6, 0) CNT at structural phase transitions. Geom- 
etry optimization of an infinite (6, 0) CNT was also performed 
by density functional theory (DFT) calculations. 



2. Methodology 

Semiempirical molecular orbital method PM3 ll2lll modified 
for one-dimensional infinite periodic structures [23] is used to 
calculate the internal energy and atomic structure of the (6,0) 
CNT under uniaxial strain. Calculations were performed using 
the MOPAC2007 code H. 
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The adequacy of the PM3 parameterization of the Hamilto- 
nian has been already demonstrated II25I1 by the calculation of 
bond lengths of the C6o fullerene with //, symmetry: the cal- 
culated values of the bond lengths coincide with the measured 
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ones [26] at the level of experimental accuracy of 10 nm at 
liquid helium temperature. The PM3 parameterization gives 
also correct bond length in graphite 0271. Comparison of dif- 
ferent types of semiempirical calculations [28] shows that the 
PM3 parameterization gives the best accuracy for structure and 
internal energy of carbon nanostructures. 

To check the convergence of calculated quantities versus the 
size of the computational cell the structure of the ground state 
of the infinite undeformed (6, 0) CNT is obtained with the use 
of computational cells comprising 72, 96 and 120 atoms cor- 
responding to three, four and five translational periods of the 
CNT, respectively. The results presented below are obtained for 
the computational cell comprising 96 atoms that allows to cal- 
culate the internal energy and bond lengths with the accuracy 
within 0.4% and 0.1%, respectively. Calculations were per- 
formed for the infinite CNT with the use of Bom-von Karman 
boundary conditions along the nanotube axis with full geome- 
try optimization (without any symmetry constraints) including 
the computational cell length. 

For the computational cell comprising 96 carbon atoms, we 
have also performed DFT calculations of the nanotube geome- 
try within package Quantum-ESPRESSO (program PWscf for 
plane wave pseudopotential method) B29I1 . Periodic boundary 
conditions along the nanotube axis with full geometry optimiza- 
tion were used, spin-polarization was not taken into account. 



PWscf program was configured to use Perdew-Zunger []30 1 
local density approximation (LDA) of exchange-correlation en- 
ergy (pseudopotential C.pz-vbc.UPF). Energy cutoff of plane 
waves was set to 408 eV. Geometry optimization was performed 
within a quasi-Newton scheme combined with the Broyden- 
Fletcher-Goldfarb-Shanno algorithm for Hessian updating (see 
also [31]). We have also checked DFT calculations by use 
of Perdew-Burke-Ernzerhof (PBE) parametrization 113 211 of 
the generalized gradient approximation (GGA) of exchange- 
correlation energy. It was found no significant difference be- 
tween the results of GGA and LDA calculations. 

In order to check the accuracy of the used DFT method we 
have calculated geometry of fullerene Ceo and compared ob- 
tained bond lengths with the experiment ll26Tl . Geometry opti- 
mization by PWscf showed certain dimerization of C-C bond 
lengths in C60 molecule. The deviation of the PWscf calcula- 
tions from the average measured values was within 1.3%, while 
for PM3 it was less than 0.2%. 



3. Ground state structure of (6, 0) nanotube 



Contrary to the armchair CNTs [6, 13], the multiplication of 
the translational period along the nanotube axis is not necessary 
for zigzag CNTs with the Kekule structure. In the work Il33ll a 
possibility of Peierls gap opening was considered in the elec- 
tron energy spectrum of the (3n, 0) CNTs in the dynamic regime 
(due to the interaction of electrons with longitudinal optical 
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Figure 1: The translational unit cell of the (6,0) nanotube, unfolded on the 
figure plane, at the ground state without deformation (a) and at relative axial 
elongation e = 10% (b); l t is the translational period along the nanotube axis 
y. Bond lengths: (a) Kekule structure: a = 0.1393 nm, b = 0.1408 nm and 
c = 0.1433 nm, d = 0.1447 nm; (b) quinoid structure: a' = 0.1449 nm and 
b' = 0.1541 nm. 



phonons at the boundary of the Brillouin zone). Since the wave- 
length of such phonons is twice as large as a translational period 
of the CNT, this leads to the doubling of the translational pe- 
riod of the system. In order to check whether the static Peierls 
distortions are realized for the (6,0) CNT, we performed PM3 
calculations with the length of computational cell equal to four 
translational periods of the CNT. These calculations as well as 
the calculations with the length of computational cell equal to 
three and five periods reveal ordinary translational period along 
nanotube axis for the (6, 0) CNT. 

Our calculations give the Kekule structure for the (6, 0) CNT 
at the ground state (s = 0), see Fig. la. There are four types of 
bonds which are arranged alternatively in the lattice. The bonds 
a and c are directed along and the bonds b and d are skewed 
with respect to the nanotube axis. The considerable difference 
of 0.005 nm approximately is found between lengths of short 
bonds a x b and long bonds c » d. 

Contrary to the PM3 method, the DFT calculations show 
only the quinoid structure of the (6, 0) CNT and do not show 
the Kekule structure for any elongation of the nanotube. How- 
ever, the DFT method seems to be less accurate in geometry 
calculations of carbon nanostructures than the PM3 method as 
being compared to experimental data (see the last paragraph of 
the methodology section). 

4. Deformational structural transition in (6, 0) nanotube 

The structure and internal energy of the (6,0) CNT are cal- 
culated with full geometry optimization for relative axial elon- 
gations of the CNT b = 0-12% (s = (I- l Q )/l , where / is the 
length of the computational cell and Iq = 1 .703 nm is the length 
of the computational cell at the ground state). The computa- 
tional cell comprising N - 96 carbon atoms with the length 
of four translational periods of CNT along the nanotube axis is 
used. 

The high kinetic barrier for formation of Stone- Wales de- 
fects excludes their spontaneous formation at room temperature 
up to the elongation s = 15% 113411 (see also consideration of 
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Figure 2: Dependence of the bond lengths a, b, c, d, a' and b' of the (6, 0) 
nanotube on the relative axial elongation e. Phase A corresponds to the Kekule 
structure, phase B corresponds to the quinoid structure. Left vertical dotted line 
indicates the lost of stability for the metastable phase B at the relative elongation 
e ss 3%. Right vertical dotted line indicates the phase transition at the critical 
relative elongation s t ~ 9%. 



these defects formation at bending deformation of CNT [35J). 
Thus, we do not consider nanotube structures with Stone-Wales 
defects. Note also that we consider CNTs at elongations which 
are well under the fracture value (e = 16% according to AMI 
calculations J36ll and e — 20% according to PM3 and density 
functional theory-based calculations H37IP . 

The following phases with different symmetries of the struc- 
ture have been found for the uniaxially deformed (6,0) CNT: 
the phase A with the Kekule structure at all considered elon- 
gations s from 0% to 12% and the phase B with the quinoid 
structure at the elongations s from 3% to 12%. The quinoid 
structure of the CNT at the great elongation is shown in Fig. lb. 
The PM3 calculated dependences of the bond lengths on the 
elongation e are presented in Fig. 2. 

The PM3 calculations of the nanotube internal energy U of 
the (6, 0) CNT at zero temperature as a function of square of 
the relative elongation is shown in Fig. 3a. Since we study the 



system at constant length and temperature, the ground state is 
determined by the minimum of the system internal energy. The 
difference of the internal energies Ua - Vb of the phases A and B 
is shown in Fig. 3b. The change of sign of the energy difference 
Ua - Ub means that the structural phase transition takes place 
at zero temperature at the critical elongation e t w 9%. Thus, at 
zero temperature, the phase A is stable at the elongation s < s t 
and the phase B (see Fig. lb) is stable at s > s t . The phase A 
at elongation s > e t and the phase B at elongation s < e t are 
metastable. 

The dependence of the bond lengths (Fig. 2) on the elonga- 
tion s demonstrates that the transition from the Kekule structure 
with four types of bonds to the quinoid structure with two types 
of bonds takes place with abrupt change of the structure of the 
CNT. Thus, first-order structural phase transition controlled by 
the nanotube elongation takes place, at zero temperature at the 
critical elongation e t m 9% where the internal energies of the 
phases coincide (see Fig. 3). At nonzero temperature, phase 
transitions in one-dimensional systems have a crossover char- 
acter (see discussion for structural phase transitions in armchair 
CNTs in [ 13] and references therein). The elastic constants (co- 
efficient of elasticity, Poisson ratio, torsion modulus) and spe- 
cific heat of the CNT also abruptly change at the phase transi- 
tion. Thus, the structural phase transition can be determined by 
peculiarities of these quantities on temperature or the elonga- 
tion. 

The dependence of the internal energy U on the length / = 
(1 +e)Iq of the computational cell is interpolated by the Hooke's 
law for the elongation range e = 0-1% 
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where v is the coefficient of elasticity of a CNT with the length 
/o = 1 .703 nm at the ground state and number of atoms N. 

Taking into account the coefficient of elasticity v defined by 
Eq. (G3, the Young's modulus of a CNT at the elongation s = 
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Figure 4: The band structure of the (6,0) nanotube with the Kekule structure (a-c), the quinoid structure (d), and the structure with all equal bond lengths (e) at 
different elongations: (a) £ = (ground state), (b) e = 2.3%, (c) £ = 8%, (d) £ = 10%, (e) £ = 0. The numbers in parts (a), (d) and (e) indicate degeneracy of 
energy bands. Only conduction bands are shown. The valence bands are the same as the conduction bands, but with opposite sign. Insets show corresponding bond 
alternation patterns. 



has the following form: 



Y = 



2nRw' 



(2) 



where R = 0.241 nm is the radius of the CNT correspond- 
ing to the elongation e = and w = 0.34 nm is the effective 
thickness of the wall 0371 - 14 ill . (In Eq. the quantity 2nRw is 
the cross-section area of a nanotube.) As a result, the value 
Y x 1.2 TPa is calculated. The obtained value of Young's 
modulus of the (6, 0) CNT with the Kekule structure is in good 
agreement with the values for this CNT with the quinoid struc- 
ture: 0.9 TPa by finite element modelling based on molecu- 
lar mechanics 113 811 . 1.42 TPa by molecular dynamics simula- 
tion using LAMMPS fH and 0.95 TPa by our PWscf DFT 
calculations. In addition, it agrees well with experimentally 
determined value 1.25 TPa for singlewalled CNTs of diame- 
ters about 1 nm ||40[ 41]. Note, that this result correlates also 



with calculations for another zigzag nanotubes: 1 . 1 TPa for the 
(10,0) CNT by PM3 calculations [37], 1.14 TPa for the (9,0) 
CNT by Hartree-Fock method Hi 1.05 TPa for the (10,0) 
CNT H, 0.94 TPa for the (10, 0) CNT H and 1.04 TPa for 
the (9, 0) CNT j44ll by density functional theory -based calcula- 
tions. 

Since the radius R depends on the elongation e, it is not con- 
venient to use Eq. (|2} for strained CNTs. Moreover, the co- 
efficient of elasticity v is the quantity which can be measured 
experimentally. Thus, we have calculated the abrupt change 
of the coefficient of elasticity at the structural phase transition: 
vb = Q.91va, where va and vg are the coefficients of elasticity 
for the phases A and B, respectively (the coefficients of elastic- 
ity are calculated for the elongation range s = 8-10%). 



same CNTs with all bonds of equal length are metallic [2C 
(Note, that tight-binding calculations, which take into consid- 
eration not only three nearest atoms but also two next-nearest 
ones [45], also give zero band gap for (6,0) zigzag CNT with 
all equal bond lengths.) Thus, the accurate calculations of the 
atomic structure of the deformed CNT are necessary to calcu- 
late the band structure. To gain a qualitative insight in the en- 
ergy spectrum of 7r-electrons in the uniaxially deformed (6, 0) 
CNT, we perform the tight-binding calculations 1E2I1 for the 
PM3 calculated atomic structure. Only interaction between 
nearest-neighbour atoms is taken into account. The hopping 
integral is used in the Harrison's form I146ll / = to(ao/acc) 2 , 
where fo = 2.6 eV is the standard value of the hopping integral 
in a CNT for the C-C bond length a = 0.142 nm 1471 148^. and 
acc is the calculated bond length. 

The band structure of the deformed (6, 0) CNT with the bond 
lengths and angles taken from the PM3 calculations is obtained 
by the diagonalization of the matrixes (lA.2t and dA.3t for s = 
0-9% and s > 9%, respectively, and shown in Fig. 4. It is 
revealed that the deformed (6, 0) CNT is semiconducting for all 
considered elongations. 

All (3n, 0) zigzag CNTs with equal bond lengths were found 
to be metallic at zero elongation B22I1 (the calculated band struc- 
ture of the (6, 0) CNT with all bonds of equal length acc = 
0.142 nm is shown in Fig. 4e). Here we have performed the 
tight-binding calculations of the band structure of the uniaxi- 
ally deformed (6, 0) CNT with equal bond lengths (that is the 
atomic structure changes only due to the changes of angles be- 
tween bonds, 

y = 2arcsin[(3(e+ l)) 1/2 /2]), 



5. Electronic structure of deformed (6, 0) CNT 

Both zigzag and armchair CNTs with the Kekule structure 
are semiconducting at the ground state with the band gap being 
proportional to the difference of the hopping (resonance) inte- 
grals corresponding to bonds of different types [20], while the 



while bond lengths remain equal. It is found that the band struc- 
ture of such CNT does not change with the elongation. 

The calculated band gap (in the center of the Brillouin zone) 
in the electron energy spectrum of the (6, 0) CNT as a function 
of the elongation s is shown in Fig. 5. For the Kekule structure 
the band gap linearly increases (from 187 meV to 311 meV) 
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Figure 5: The band gap E g in the ^-electron energy spectrum of the uniaxially 
deformed (6, 0) CNT versus its axial elongation e. (The solid line is a guide to 
the eye). 

with the elongation for the elongation range s = 0-2.3% and 
linearly decreases (from 311 meV to 88 meV) for the elonga- 
tion range e = 2.3-9%. The bond lengths linearly increase with 
the elongation. The phase transition from the narrow band gap 
(« 0.1 eV) to the moderate band gap (« 0.5 eV) semiconductor 
occurs in the electron subsystem of the CNT at the elongation 
s t as 9% corresponding to the structural phase transition. Un- 
like the Kekule structure, the band gap of the quinoid structure 
increases (from 513 meV to 730 meV) with the elongation for 
the considered elongation range s = 9-12%. 

6. Strain sensor based on zigzag CNT 

At present, the standard methods of the nanometer-scale de- 
tection of strain are based on micrometer size silicon piezoresis- 
tors or MOS transistors with piezoresistive channels [49]. Re- 
cently, a strain sensor based on commensurate-incommensurate 
phase transition in the double-walled CNTs was proposed |01 ■ 
Since the resistance of an undoped (intrinsic) semiconductor 
depends exponentially on the band gap, the predicted effect of 
the abrupt band gap change at the structural phase transition 
controlled by the axial elongation of the (6, 0) CNT can be de- 
tected by measurements of the CNT resistance. Thus, a strain 
sensor based on measurements of the resistance of the uniax- 
ially stretched (6, 0) CNT can be proposed. The preliminary 
stretched (6, 0) CNT with the initial elongation £j n i can be at- 
tached to electrodes which are positioned on controlled object 
(see Fig. 6). If this object is compressed or extended, the elon- 
gation of the CNT is s — s m i + e Q bj, where e bj is the object 
elongation. At the object elongation e Q bj = fit _ £ ini, which cor- 
responds to the phase transition in the CNT, the abrupt change 
of the resistance can be detected. The proposed strain sensor 
can be used for detection both compression and extension of 
the object for the initial elongations of the CNT s m i > s t and 
£ini < s t , respectively. Since the value of a typical diameter of 
a single-walled CNT is within several nm, a minimal size of a 
strain sensor is limited mainly by a minimal possible width of 
the electrodes. 

The (6, 0) CNT is considered here as an example of (3n, 0) 
zigzag CNTs. We believe that the phase transition from a nar- 
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Object to be controlled (insulator) 

Figure 6: The nanometer-scale strain sensor based on the stretched (6, 0) CNT. 
If initial nanotube elongation is less than critical one £ m \ < e t , then the strain 
sensor operates for registering an extension of an object (e bj > 0), and if > 
St, it operates for registering a compression of an object (e (,j < 0). 

row gap to a moderate gap semiconductor is possible also at 
elongation of other (3«,0) CNTs. Thus, the proposed strain 
sensor can be based also on (3«, 0) CNTs as well as on (n, n) 
armchair CNTs, where a narrow gap semiconductor to metal 

J 1 

transition is predicted at elongation 111 3D . A set of (3n,0) or 
(«, n) with different chiral indices (and thus with different criti- 
cal elongations e t ) or with identical chiral indices but with dif- 
ferent initial elongations can be attached to the same pair of 
the electrodes. For such a scheme of the strain sensor, the tran- 
sitions take place for different CNTs at different elongations of 
the controlled object, therefore the sensor can be used not only 
to detect the presence of strain in the controlled object but also 
to measure the value of this strain. 

7. Conclusions 

The first order deformational structural phase transition at 
zero temperature is revealed at the critical uniaxial elongation 
e t « 9% of the (6, 0) carbon nanotube by semiempirical PM3 
molecular orbital calculations. Namely, it is found that the (6, 0) 
nanotube has the Kekule structure with four types of bonds and 
the quinoid structure with two types of bonds for the elonga- 
tions less and greater than the critical value e t > respectively. 
Geometry optimization by the semiempirical PM3 method is 
found to be better than the DFT results as compared to exper- 
imental data on the C6o molecule. Tight-binding band struc- 
ture calculations show that this structural phase transition is 
followed by the transition in the electron subsystem from the 
narrow gap (w 0.1 eV) to the moderate gap (« 0.5 eV) semi- 
conductor. A strain sensor of nanometer size based on mea- 
surements of resistance of the uniaxially stretched nanotube is 
proposed. 
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Appendix A. Tight-binding model details 

A.l. Kekule structure 

For the Kekule structure with four types of bond (a, b, c and 
d) the translational unit cell (Fig. la) has three times more in- 
equivalent atoms than the quinoid structure with two types of 
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Figure A. I : The two-dimensional primitive unit cell (shown by dotted lines) of 
the zigzag carbon nanotube with the Kekule structure. The coordinate system 
x, y is used for calculations of the Hamiltonian matrix JA.2t . 




Figure A.2: The two-dimensional primitive unit cell (shown by dotted lines) of 
the zigzag carbon nanotube with the quinoid structure. The coordinate system 
x, y is used for calculations of the Hamiltonian matrix IA.31 . 



bonds a' and V (Fig. lb). Therefore, for the Kekule structure 
the tight-binding Hamiltonian matrix has dimensionality 6x6 
which is three times greater than the one for the quinoid struc- 
ture. 

In the case of the Kekule structure with four different 
bond lengths (Fig. la) the two-dimensional primitive unit cell 
(Fig. A.l) contains six carbon atoms (three a and three /3 
atoms). 

By employing the tight-binding approach for the Kekule 
structure, a Bloch wave function can be written as HUH: 

¥(r) = -jL V exp(ikR c ) ^(r - R ( ), (A. 1) 

where N t is the number of translational unit cells, R^ is the po- 
sition of the £-th kind of atoms = a\,a%, a^jSi./fej/fe), and 
iff z (r) is the normalized p z atomic orbital wave function of an 
isolated carbon atom. 

In the tight-binding approximation, the calculations of a band 
structure of a CNT with the Kekule structure are reduced to 
diagonalization of the six order matrix. The Hamiltonian matrix 
is obtained as (see also fjoll ): 



(A.2) 



where t a , %, t c , tj are the hopping integrals corresponding to the 
bonds a, b, c, d, respectively; k = k x e x + k v e v is the electron 
wave-vector; e c , e v are the unit vectors of the coordinate system; 
ri = aifii = b x e x -b y e y , r 2 = anJ3 2 = -d x e x -d y e y , r 3 = ary8 3 = 
c y e y = ce v , r 4 = a 2y 8 2 = UiPi - d x e x - d y e y , r 5 = ff 2 /? 3 = 
-b x e x -b v e v , = affi = a y e y - ae y are the vectors connecting 
nearest c,- and atoms; a x = 0, a y = a, c x = 0, c y = c, 
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-bcosS, d x = c/siny, d y = -of cosy are the 
projections of the bonds on the x and y axes of the Cartesian 
coordinate system; y and 6 are the angles between the bonds a 
and d, and the bonds b and c, respectively. 



The periodic boundary condition is 
Ch-k = 2jrq; q = 1, . . 
nq 



k x 



2nq 

"c7 



.,2n/3, 



n , n 



2{2d x + b x ) 

where Ch = 4(2*4 + b x ) is the magnitude of the chiral vector 
Ch = na\ + 0a 2 = (n, 0) (for n — 6 the value of Ch is a circum- 
ference length of the (6, 0) CNT; ai and a 2 are the unit vectors 
of graphene), and l t -a + c + 2d y is the period of the (6, 0) CNT 
with the Kekule structure along the nanotube axis. 

A.2. Quinoid structure 

In the case of a CNT with the quinoid structure with two dif- 
ferent bond lengths (a' and V as shown in Fig. lb) the primitive 
unit cell contains two atoms (denoted as a and /3 in Fig. A.2). 
In this section we will denote bonds a' and V as a and b, re- 
spectively. 

The Bloch wave function for the quinoid structure can be 
written in form of Eq. (IA.U . where R( stands for the positions 
R a and Rg of the a and /3 carbon atoms. 

Th e tight- binding Hamiltonian matrix in this case has the 
form|2a|22| 

H% 

Ha/3 



H 



(A3) 



,ikr, 



+ f fl e ikr3 and W p = t h e 



ikr, 



where H a/3 = t b e' kr < + t a & 
t a e- [kt2 +f fl e- ikr -'. 

As it follows from Fig. A.2 the vectors ri, r 2 , r3 are defined 
in the chosen reference system as: 

ri = be y , r 2 = -a x e x - a y t y , r 3 = a x e x - a y e y , 

where a x = asiny, a y = -a cosy are the projections of the 
bond a on the x and y axes; the wave-vector is k = k x e x + k y e v . 
The periodic boundary condition is 



Ch-k = 2nq; q = 1, 
2nq nq 



,2n, 



, 7i ^ ^ n 

Ch 6a x ' l t l t ' 



where Ch = 2na x — I2a x is the magnitude of the chiral vector, 
and l t = 2(a y +b) is the period of the (6, 0) CNT with the quinoid 
structure along the nanotube axis. 

To obtain the energy spectrum the following secular equation 
should be solved 



det 



-E H 



0. 



By expanding the above determinant, we obtain 

E 2 - t 2 b + At 2 a cos 2 (nq/6) + 4t a t b cos(k y l t /2) cos(nq/6). 

Thus, in the tight-binding approximation the dependence of 
7r-electron energy on the wave vector along the nanotube axis 
E m (k) (band structure) for a («, 0) CNT with the quinoid struc- 
ture is [20]: 



EJk) = + [t 2 b + At 2 a cos 2 (nq/n) 

+ At a % COS(7T<7/«) cos(W t /2)] 



1/2 



(A.4) 



where t a and tf, are the hopping integrals for the bonds a and b, 
respectively, q = 1,2, . . . ,2n is an integer, and -n/l t ^ k y ^ 
n/l t is the wavenumber of the 7r-electron along the nanotube 
axis y. 

For the CNT with the quinoid structure there are An energy 
bands in the electronic spectrum of the CNT, An- A of which are 
twofold degenerate, according to equivalence of two possible 
directions of the electron quasi-momentum. Only bands with 
indexes q equal to n and 2n are not degenerate. The signs plus 
and minus in the formula ( IA.41 > correspond to the conduction 
and valence bands, respectively, which are symmetric relative 
to the Fermi level E — 0. For the (n, 0) CNT with even «, bands 
with index q = n/2 are dispersionless (the electron energy does 
not depend on the quasi-momentum) (see also [22]). The for- 
mula ( IA.4b for the electron energy bands shows that the elec- 
tron energy in the (n, 0) CNT can be equal to zero only when 
n is multiple of 3 and t a = % (i.e., all equal C-C bonds). Only 
in this case the electric conductivity of the CNT has a metallic 
character, in all other cases there is a gap in the electron energy 
spectrum. 
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